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ABSTRACT 

Properties  of  two  alternative  procedures  to  the  Jack- 
knife  Point  and  Confidence  Interval  Estimation  Procedure  of 
Gaver  and  Chu  have  been  studied.  They  are  called  the  Log- 
Normal  Likelihood  Procedure  (LNLJ)  and  the  Moment  Procedure 
(MP) .   These  two  procedures  were  investigated  and  compared 
with  the  Jackknife  Point  and  Confidence  Interval  Availabil- 
ity Estimation  Procedure.   Numerical  results  from  simula- 
tions are  presented  in  this  report. 
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I-   INTRODUCTION 

A.   OVERVIEW 

Availability  is  the  measure  of  equipment  effectiveness 
that  relates  reliability  and  maintainability  to  operational 
readiness.   In  some  cases,  availability  and  operational 
readiness  have  been  considered  the  same.   In  general,  these 
are  all  requirements  which  must  be  considered  during  the 
design  and  operation  of  a  system,  and,  thus,  must  be  quanti- 
tatively evaluated. 

There  are  a  variety  Df  ways  of  expressing  availability. 
In  general,  availability  relates  up  time  (reliability- 
related)  to  down  time  (maintainability-  related) ,  and  it  may 
be  defined  as  the  ratio  between  the  time  the  system  is  capa- 
ble of  performing  its  mission  to  the  total  time  the  system 
is  in  operational  demand.   Alternatively,  it  is  a  measure  of 
the  probability  that  a  single  system  is  "up"  or  available  at 
a  possibly  "random"  point  in  time  when  its  services  are 
needed. 

The  three  expressions  of  availability  of  greatest  con- 
cern are  --  (1)  inherent  availability  (Ai)  ,  (2)  operational 


availability  (Ao) ,  and  (3)  achieved  availability  (Aa) .   See 
Rise  and  B  jorklund[  Ref .  1  ].   In  this  study,  we  will  be  con- 
cerned with  only  the  inherent  availability.   Often,  inherent 
availability,  which  is  a  hardware  oriented  measure,  is  spe- 
cified and  required  within  a  maintainability  contract 
requirement. 

A  matter  of  primary  concern  is  that  of  estimatinq  avail- 
ability from  observations  on  system  up  times  and  down  times. 
The  estimating  process  should  furnish  both  a  point,  or  sin- 
gle number,  estimate  and  also  some  measure  of  the  stability 
of  the  estimate  such  as  a  standard  error,  or  confidence  lim- 
its.The  latter  problem  is  more  difficult  than  the  former. 

In  the  paper  by  Gaver  and  ChuTRef.  2],  it  has  been 

demonstrated  that  the  jackknife  technique  can  be  useful  for 

confidence  interval  estimation  of  system  availability .For 

instance,  let 

2(0) 

A  = (1.1) 

3 (U) *B(D) 

be  the   long   run   system  availability    of   a    single   unit    system 

that  changes   from    up   to   down  states    in   accordance    with   a 

two-state   semi-Markov    (or   more    general)     process   for    which   up 

[down]    state   expected   duration    is    3(0")  [3(D)].      Then    it    is 


possible  to  successfully  jackknife  the  naive  point  estimator 

that  involves  replacing  expectations  by  the  corresponding 

sample  means, 

0 
"J-  0.2) 

5  ♦  D 
after  initial  logistic  transformation.  There  are,  however, 
particular  parametric  families  of  distributions,  popular  as 
models  for  summarizing  up  and  down  time  data,  for  which  the 
most  statistically  efficient  estimators  of  E  (U)  and  2(D)  are 
not  the  simple  sample  arithmetic  means. In  particular,  this 
is  true  for, 

(1) .The  log-normal  distribution   with  unknown  (to  be 

estimated)  parameters;  often  this  is  used  as  a  model 
for  down  or  repair  times.   In  particular,  the  log 
normal  distribution  has  enjoyed  considerable  popu- 
larity for  representing  electronic  system  repair 
times. See  Kline  and  Almog[Ref.  3]. 
(2). the  general  gamma  distribution  with  unknown  shape 

and  scale  parameters,  and 
(3). the  general  Weibull  distribution  with  unknown  shape 
and  scale  parameters. 
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B.   PURPOSE  AND  APPROACH 
1 .   Objective s 

Availability  estimation  by  use  of  the  jackknife  has 
already  been  studied  by  Saver  and  Chu.This  thesis  presents 
two  procedures  alternative  to  the  above.   These  may  be 
called  the  Log-normal  Likelihood  jackknife  (LNLJ)  procedure 
and  the  Moments  procedure  (MP).   In  brief,  the  jackknife 
method  has  the  capacity  to  reduce  the  statistical  bias  of 
estimates  of  such  quantities  as  system  availability,  and 
also,  and  more  importantly,  to  furnish  confidence  limits 
that  behave  in  a  satisfactory  manner  (economically  enclose 
the  true  availability)  despite  the  fact  that  underlying  dis- 
tributions are  unknown.   See  Miller[Ref.  4]  for  a  review  of 
much  of  the  literature  on  jackknifing. 

The  present  estimation  procedure  makes  use  of  the 
specific  assumption  of  a  log  normal  model  for  repair  times 
to  compute  the  maximum  likelihood  estimate  (mle)  of  avail- 
ability.  Properties  of  the  procedure  have  been  studied  by 
Monte  Carlo  simulation.   A  number  of  such  simulation  results 
are  presented  in  this  thesis.   Also  presented  are  compari- 
sons with  the  direct  method  of  Gaver  and  Chu,  and  with  a 
computationally  simple  method  based  on  moment  estimation. 
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Attention  also  has  been  paid  to  study  of  the  sensitivity  of 
the  LNLJ  procedure  to  specification  error.   This  refers  to 
the  error  introduced  by  assumption  of  a  log-normal  distribu- 
tional model  and  the  corresponding  maximum  likelihood  esti- 
mates when,  in  fact,  another  distribution  governs  the  obser- 
vations.  It  is  found  that  LNLJ  method,  which  is 
theoretically  most  efficient  in  large  samples  if  the  model 
(likelihood)  is  correct,  can  sometimes  produce  consistently 
biased  results  when  the  log  normal  model  does  not  apply. 
Since  it  is  nearly  impossible  to  assure  log  normality  from 
small  samples  of  data,  this  finding  suggests  that  caution  is 
in  order. 

2.   L03  normal  Likelihood  Procedure  For  a  Single  On it 
(a). Data  is  given  as  follows, 

u1,u2, ,un  (up  times) 

d  ,d  , ,dn  (down  times) 

x1=lnd1  ,x2=lnd2  , ,xn  =  lndn 

(b). Assume  down  times  are  log-normal. Then, 

E[D]«  exp  (  nd+  ^/2•   ad*)  (1.2.1) 

where 

u  =  E[lnD]  (1.2.2) 

a 

a2  =  7ar[lnD]  (1.2.3) 

d 
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Furthermore,  the  maximum  likelihood  estimates 
(mle)  of  y,  and  a2  are 
^d  =  x  =  1/n  2  ln(d.)  (1.2.4) 

a2   =  1/n  I   (x .-x)  2  ,  (1. 2.5) 

d  1- 

so   the  mle   of    E[D]    is,    by    invariance, 

E[D]    =  exp(ud*    1/2*  7|)  (1.2.6) 

it  is  essentially  the  latter  expression  that 
will  be  used  in  the  estimate  for  availability  to 
replace  the  simple  first  moment  estimate  E(D)=d. 
In  fact,  v  _     has  been  replaced  by  its  unbiased 

d 

version    (n   replaced   by    n-1). 
(c).If  the  above   holds,    it   is  advantageous   to   trans- 
form   first    (See    Hosteller    and   Tukey[ Ref .    5])    the 
estimated   availability 
S(O-) 


A  = 


B(D)  +E  (D) 
the  log-logistic  transform  of    availability   is 

ln(l   /   1-A)    =    lnE*(u)-(ul +1/2*? /)  (1.2.7) 

a  d 

Jackknifing  will  be   carried   out    using    the  sta- 
tistic 
L     =    In  (\   /   1-A)    =ln(u)-  (x+1/2«sx2)  (1.2.9) 
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(d) .Recompute  L   repeatedly,  leaving  out  succes- 

LiN 


sively  the   sample   pairs    (u    ,d    ),    (u    ,d    ),    , 

11  2       2 

{u    ,d   ) 

3-1  fl 

u  =    1/n-L  (£  a  .+£ui   }  (1.2.9) 

J  i  =  l  »«j+» 

x.jVn-1»{^  x.+  J  x,     }  (1.2.10) 


s   2  =  1/n-2»f  ^(xi  -   X    )2  +    ^  (x^-    x   )2    }  (1.1.11) 

1^   j    InM-^j*    V2.sxaj  (1.2.12) 


where    j=1  ,2, , n 

(e). Compute   the  pseudovalaes   as   follows: 

PLN,T  nLLN-(r--1)I'LN,-j  J'1'2 C      <K2-13» 

Recall  that  L   =L    is  the  result  of  the  com- 

LN   LN,all 

puting  the  quantity  to  be  jackknif edr  leaving 
out  none  of  data. 
(f). Compute  the  mean  and  variance  of  -he  pseudova- 
lues, 


LN 


J-.  n 


S  2=  1/n-1«  T   (P   -  P    2  (1.2.15 

LN  *-    LN.j   LN 

(g) .The    jackknifed    point  estimate    of   the    availabil- 
ity  is   now, 
A=    exp(PLN)     /    1+exp(PLN)  (1.2.16) 

(h) . Symmetric   two   sided   confidence    limits   at    confi- 
dence   level    (!-*)%   are   derived   as    follows; 

CL      -    P       -    t  Cn-1)     5QST(1/n.S    2)  (1.2.17) 

LN  LN  l-ct/2  LN 
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C%=    *m+   tl-a/2(n-1'    s2ET(1/r».SJ)  (1.2.18) 

where   t.      /0(n-1)    is   the    (1-a/2)100S    quantile   of 

l-a/2 

Student's-t  with  n- 1    degrees   of    freedom.      Then, 

exptCL^)                            axp(CULN) 
—    <   A    < '  (1.2.19) 


1    +    exp(CL      )  1    ♦    exp(CU      ) 

LN  LN 

with   confidence   approximately    (1-ot  )  100%.       Note 
that    the    confidence    limits  are   nearly   symmetric 
around  In  (E[0  ]/E[D])  ,    and   not   around    A. 
(i) .One   sided   confidence   limits   at    confidence   level 
(1- a  )100%   are    derived 

CLLN=    *LN~    Va    (n"1)     SQRT(1/n.Sj)  (1.2.20) 

C0LN=    ?LN  +    *l-a    <n"1)     SQRT(1/n.SL2)  (1.2.21) 

So,    one-sided   tipper   confidence   limit    is 

expfCU^) 

A   <     (1.2.22) 

1    *    exptCO^) 

and   lower    confidence   limit   is 

expfCL^, 
A    <    (1.2.23) 

1    *    exp  (CL^) 

Note   that    if   up  times  are   exponentially   distributed 
then   u    is   actually  the   mle    of  E[U],    so   under   the    model 
assumptions   the    (transformed)    mle   of    availability    is    being 
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jackknif ed ;    this    procedure    has    been   validated   by    Reeds 
[Ref.    6]   for    large   samples.      Simulation   studies   are 
nevertheless  essential   for    investigating    properties    of   the 
jackknife  or   alternative   techniques   for   small  samples,    and 
for   investigating   sensitivity  to   model   assumption. 
3.       Moment   Procedure   For   a    Sin^il   Unit 
(a). Data    given   as    follows 


VU2' 


,a. 


(up  times) 
(down  times) 


71=lnu1,y2=lnu2,....#yn=lnun 

x1  =lndx  ,  x2  =lnd2  , . . .  .  ,y^  sladn 
(b). Assume   up    and    down    times    are    log-normal. Then, 

1. 3.  1) 
1.  3.2) 


E[U]=    exp(uu*    1/2. au2) 


E[D]=    exp  (ud+    1/2»ad2) 
where, 

uu=    E[lnO] 


u    =    E[lnD] 

a 

a2=    7ar[ln0] 

u 

a^=    7ar[inD] 


Furthermore,    the  mle   of  p     ,   az,  y     ,    and  a 

u  u  d 

...    /»< 


y 


~u  =  1/n#  X  (7i"y} 


td=  x   =1/n.  £  r. 


[1. 

3. 

3) 

1, 

3. 

») 

[1, 

3. 

5) 

;i« 

3. 

6) 

2 

d 

are 

[1, 

3. 

7) 

1, 

3. 

3) 

1. 

3, 

9) 
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:  =    1/n«    £  (x.  -x)  2  (1.3.  10) 


2   ^ 
d 


so   the  ile   of   E[U]    and    E[D]   are,    by    invariance, 
E[U]=    exp(yu+    1/2*^2)  (1.3.11) 

E[D]=    exp(u*    1/2*  o*)  (1.3.12) 

a  a 

(c)  If  above  holds,  the  logistic  transform  of  avail- 
ability is, 

In  (A  /  (1-1))  =  (T^+1/2  "u2)-(ud+1/2  ad2)  (1.3.13) 

estimation   will    be   carried   oat    using   the   statis- 
tic 
P»   In  (A   /    (1-A)  )  =  (y  +  1/2»sy2)-(x  +  l/2.sx2)  (1.3.14) 

(d)  Simply  the   point  estimate    of   availability   is 
now, 

A=   exp(P)    /    H-exp(P)  (1.3.15) 

(e)  Variance    of  the    statistic   P   can   be   calculated   as 
follows , 

Var  (P)  =  (SE)  2=    i/n«sy2    ♦    1/4«Var[ sy2)    ♦    1/n«sx2    ♦ 

1/4«Var[s    2]  (1.3.16) 

where, 

7ar[sy2]=    i/n   [M4y-(Sy2)2]  (1.3.17) 

7ar[sx2]=    1/n   [M4x-(sx2)2]  (1.3.13) 

(f)  Symmetric   two   sided   confidence   limits   at    confi- 
dence  level    (1-a)%   are   derived   as   follows, 

L=    P    -    Z1_c/2(SE)2  (1.3.19) 
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0=    P    *    *L-a/2  (SE)2  (1.3.20) 

where    z1_a/2is    the    (1-  a/2)  100%    quantile   of   the 

standard   normal    distribution.    Then, 

(exp(L)     /    Uexp(L)}    <   k    <    (exp(O)    /    H-exp(U)} 

(1.3.21) 
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II.   SIMULATION  PROCEDURE 

A  simulation  procedure  has  been  used  to  compare  the 
operating  characteristics  of  LNLJ  methodology  introduced 
here  with  other  approaches.   Specifically,  simulation  has 
been  used  to  compute 

(a)  the  actual  coverage  of  the  true  availability  figure.  A, 
by  the  confidence  intervals  given  by  the  procedures 
under  study,  when  the  nominal  coverage  is  (1-  a  ) 100^, 

(b)  measures  of  confidence  interval  sxpected  width  and  vari- 
ance width, 

(c)  estimate  of  the  expected  point  availability  estimated  by 
the  procedure  under  study. 

The  simulation  program  was  written  in  FORTRAN  IV,  and 
the  simulation  have  been  carried  out  on  the  IBM  3033  at  the 
Naval  Postgraduate  School.   A  more  detailed  description  of 
the  program  and  its  major  subroutines  is  given  in  Appendix 
A. 

An  outline  of  the  simulation  procedure  now  follows; 

(A)  A  simulated  sample  of  up-rime  durations,  (ui#  i=i, 
2,...,  n)  and  down-time  durations,  (di ,  i=1,  2,...,n)  were 
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times.      The    Naval    Postgraduate   School   LLRANDOM    package   was 
used,    along   with   the   International   Mathematical   and    Statis- 
tical  Library    (IMSL)    random   number    generator.      For    all   cases 
considered     (exceptions   noted) ,    expected   up    and   down   times 
were   E[U]=73.036    and   E[  D  ]=3 .  8U4  231    for   the    population   sam- 
pled,   so   long-run    availability    A=0.95.      Sample    sizes    of   n=15 
and   n=25    were   simulated,    and  provided  the   bases   for    point 
estimators    and    for  confidence   intervals. 

(B)     Here  are   the   different    distributional   situations 
that   were   sampled.      A   total    of    1000    replications   was    used  to 
evaluate  each   procedure  in    each   distributional    situation. 

(1)  Dp-times    come    independently    from   the   iid   Expon(\) 
distributions.      Down-times   come    independently   from 
the   log-normal  (y, a2)     distribution. 

(2)  Up-times    are   iid   Expon(\).      Down    times    are   iid 
Expon  (u)  . 

(3)  Op-times    are   iid   Expon(\).       Down-times    are    iid 
Gamma  (9, )c),         being    a   scale   and   k   a   shape    parameter, 
so    arranged   that    E[  D  ]=  f*~l    and   7ar[D]=(/T^  )~2.      A 
Gamma   with   k>1,    e.g.    k=2,    represents   data    that    is 
more   tightly    grouped  around   Lzs   mean   than    is    true   of 
exponentially   distributed   data;    it    roughly   resembles 
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data  that   is   log-normal.      A   gamma    with   k<1,    e.g. 
k=1/2,    represents    data    of   relatively  extreme   disper- 
sion and    positive    skewness   as   compared    to    exponen- 
tial  data   having   a    very    long    right    tail,    compensated 
for    by   a    high   density   near   zero.       A   gamma    distribu- 
tion  with   integer   k     (k= 1 ,2,  3, . . . . )    is  often   called 
Erlang;    realizations  are  easily   simulated    by    summing 
k   independent   Exponential   realizations. 

(4)  Op-times   are   iid   Expon  (A).      Down-times   are    iid 
long-tailed   h   distribution: 

d=E[D]«w«exp(h«w)»  (1-h)  2    ,         0<h<1, 
where  w   comes   from   Expon(1).      This   long-tailed   h 
distribution    possesses    no  closed-form   representa- 
tion.     However,    it    can    be   shown   to    have   characteris- 
tics  similar   to  those   of   the   Exponential    for    small 
values   of   h,    but    to    have  a   systematically    longer 
right   tail  than  the    exponential. 

(5)  Up   times    are   iid   long-tailed   h  distribution-      Down 
times   are   iid   Expon  (ji) . 

(6)  Op    times    are    iid   Expon (a).      Down    times    are    iid   Wei- 
bull  (6, k),         being    a  scale    parameter   and    k    a    shape 
parameter,    so   arranged    that   S[D]=9l/k*  r(1/k    +    1)    and 
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Var[D]=e2/k*[  r  (2/k   ♦   1)-  r2OA   *    1)  ].      A    weibull 
with   k>1,    e.g.    k=2,    represents  data   that    is   more 
tightly  grouped   around    its    mean   than   is   true   of 
exponentially   distributed   data;    it    quantitatively 
resembles   data  that    is    log-normal. A   Weibull    with 
k<1,   e.g.    k=1/2,    represents   data   of   relatively 
extreme   dispersion    and    positive  skewness   as   compared 
to    exponential  data    a   very   long   right  tail,    compen- 
sated  for    by   a   high    density    near    zero. 
(7)     Op    times    are   iid   Sxpon(M  .      Down    times   are   iid 

long-tailed   log-normal   h   distribution.      See    Appendix 
B. 
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III.   ANALYSIS 

The  methods  described  were  simulated  for  various  distri- 
butional assumptions  which  were  defined  in  the  Chapter  2. 
Simulation  results  for  each  method  and  set  of  distributional 
assumptions  are  shown  in  Table  1  and  Table  2. 

In  general,  HP  has  very  high  coverage  factor  and  satis- 
factory point  availability  estimation.  However,  both  average 
and  variance  of  confidence  length  are  higher  than  JK  and 
LNLJ  procedures. In  addition,  no  consistent  superiority  of 
LNLJ  procedure  over  JK  procedure  has  been  noted. 

Especially  under  the  Exponential  up  and  Garana  (k=1/2) 
down,  and  Exponential  up  and  We ibull  (k=1/2)  down  times 
assumptions,  a  consistently  low  coverage  factor  and  point 
availability  estimation,  and  a  consistently  high  average 
confidence  length  and  variation  have  been  obtained  from  LNLJ 
procedure.   In  these  cases,  distributional  characteristics 
of  down  times  after  log-transform  provide  highly  left  skewed 
shape  and  large  variance,  and,  also  these  features  reflect 
in  the  pseudovalues  which,  in  turn,  causes  left  skewed  shape 
rather  than  a  symmetric  shape  for  the  pseudovaiue 
distribution. 
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In  order  to  compensate  for  this  situation,  the  Biweight 
procedure,  see(Ref.  7  ],  and  the  Winsorizing  procedure,  see 
[Ref.  8],  were  implemented  directly  on  the  pseudovalues. 
Results  are  shown  in  Tables  3,  4,  5,  6,  7,  8.   These  proce- 
dures are  robust/  resistant  methods  for  establishing  confi- 
dence limits  on  means;  they  tend  to  down-weight  extreme 
observations  that  appears  as  outliers. 

For  small  sample  size,  (i.e.  n=15)  Winsorizing  at  g=2 
and  g=3  levels  provide  an  adequate  coverage.  3ut,  the  point 
availability  estimate  and  average  confidence  length  remain 
inadequate.   In  general,  the  Biweight  statistical  procedure 
caused  reduction  of  the  coverage  factor,  and  it  is  not  yet 
considered  effective  for  compensation. 

Next,  a  grouped  jackknife  procedure  was  applied  to  these 
cases. Simply,  up  and  down  times  data  were  grouped  two  by 
two,  and  average  of  these  groups  have  been  taken.  Obviously, 
this  process  reduced  the  sample  sizes  to  one  half  and  caused 
loss  of  degrees  of  freedom. However,  smaller  variation  and 
more  stable  results  were  anticipated. Simulation  results  are 
presented  in  Tables  9  and  10.   Very  good  coverage  factors 
and  point  availability  estimates  were  obtained,  but,  average 
confidence  length  remained  at  a  rather  high  level.   See 
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Appendix  C  for  histograms  of  both  cases  for  up  and  down 
times,  and  log- transform  of  down  times  and  pseudovalues 
after  grouped  jackknife  procedure  applied. 

Other  critical  simulation  results  were  obtained  from 
exponential  up  and  long-tailed  log-normal  h  down  times 
case. For  both  JK  and  LNLJ  procedures,  above  statistical 
procedures  were  applied,  and  results  were  shown  in  Tables 
11,  12,  13,  14,  15,  16. In  either  case,  after  implementation 
of  these  statistical  representations,  the  coverage  factor 
decreased.   Also,  see  Appendix  C  for  histograms  of  this 
case. 
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Table    1:      Simulation   Results   for   Distributional    Assumptions 

Sample   si2e=15 

Underlying  Average  Variance  Point 

Distribution             Coverage  Width  Width        Availability 

A.  Exponential    up      JK      0.9380  0.0838  0.0022  0.9458 
Log-normal   down    LN      0.9443  0.083^  0.0016  0.9459 

MP      0.9698  0.0959  0.0026  0.9556 

B.  Exponential    up      JK      0.9470  0.0850  0.0013  0.9471 
Exponential   downLN      0.9417  0.1378  0.0123  0.9312 

MP      0.9868  0.1765  0.0196  0.9421 

C.  Exponential    up      JK      0.9423  0.0706  0.0007  0.9465 
Gamma   down               LN      0.9419  0.0764  0.0010  0.9437 

k=2                          MP      0.9830  0.1001  0.0031  0.9538 

D.  Ditto                     JK      0.9450  0.1131  0.0039  0.9442 
k=1/2                     LN       0.9161  0.4574  0.0932  0.8064 

MP      0.9873  0.5013  0.0970  0.83  11 

E.  Exponential    up      JK      0.9256  0.1171  0.0097  0.9437 
Long-tailed   h        LN      0.9483  0.1750  0.0197  0.9287 

h  =  0.2                     MP       0.9780  0.1951  0.0244  0.9412 

F.  Ditto                     JK      0.3753  0.1663  0.0343  0.9387 
h=0.4                       LN       0.9406  0.2144  0.0342  0.9314 

MP       0.9513  0.2092  0.0321  0.9431 
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S.    Long-tailed   h    upJK      0.9245        0.1077      0.0024  0.9428 

Exponential   downLN      0.9058        0.1650      0.0150  0.9256 

h=0.2  MP      0.9803         0.1973      0.0212  0.9381 


H. 


Ditto 


h=0.4 


JK      0.8837         0.1470      0.0048 


LN      0.3552         0.2100      0.0179 


MP       0.9437         0.2377      0.0244 


I.    Exponential    up      JK      0.9476        0.0641      0.0005 
Weibull    down  LN      0.9449        0.0672      0.0007 


k  =  2 


MP      0.9838        0.0917      0.0023 


0.9351 
0.9163 
0.9284 

0. 9483 
0.9463 
0.9501 


J. 


Ditto 


k=1/2 


JK      0.9142         0.1764      0.0239 


LN      0.9336         0.5752      0.0921 


MP      0.9751         0.5319      0.0966 


0. 9404 
0.76  49 
0. 8464 


K.  Exponential  up  JK  0.8968  0.1385  0.0220 
Long-tailed  log  LN  0.9119  0.1217  0.0091 
normal   h    down        MP      0.9183         0.1113      0.0058 


L. 


h=0.2 


Ditto 


h  =  0.4 


JK      0.8794         0.1578       0.0305 


LN       0.9009         0.1394      0.0136 


MP       0.9001         0.1210      0.0079 


0.9416 
0.9483 
0.9576 

0.9394 
0.9483 
0.9575 
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Table   2:      Simulation   Results   for   Distributional   Assumptions 

Sample   size=25 

Underlying  Average  Variance  Point 

Distribution  Coverage      Width  Width        Availability 

A.    Exponential   up      JK      0.9436  0.0594  0.0007  0.9470 

Log-normal   down   LN      0.9U73  0.0589  0.0004  0.9471 

MP      0.9740  0.0713  0.0008  0.9571 

3.    Exponential    up      JK      0.9473  0.0597  0.0004  0.9471 

Exponential    downLN      0.9228  0.0978  0.0052  0.9322 

MP       0.9945  0.1302  0.0086  0.9450 

C.  Exponential    up      JK      0.94  75  0.0506  0.0002  0.9474 
Gamma    down                LN      0.9414  0.0547  0.0003  0.9447 

k=2                          MP      0.9895  0.0751  0.0010  0.9552 

D.  Ditto                     JK      0.9464  0.0766  0.0009  0.9463 
k=1/2                     LN      0.8187  0.4008  0.0790  0.8168 

MP      0.9884  0.4506  0.0841  0.8435 

E.  Exponential    up      JK      0.9313  0.0834  0.0041  0.9456 
Long-tailed    h         LN      0.9464  0.1213  0.0  077  0.9323 

h=0.2                      MP      0.9867  0.1426  0.0108  0.9450 

F.  Ditto                     JK      0.3734  0.1261  0.0212  0.9421 
h=0.4                     LN       0.9475  0.1464  0.0132  0.9366 

MP      0.9593  0.1510  0.0139  0.9481 
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G.    Long-tailed   h    upJK      0.9297        0.0759      0.0007 

Exponential   downLN      0.8901         0.1167      0.0066 

h=0. 2  MP      0.9834         0.1447      0.0099 


H. 


Ditto 


h=0.  4 


JK       0.8802         0.1053      0.0017 


LN      0.8248         0.1521       0.0073 


MP      0.9575         0.1771       0.0116 


I.    Exponential    up      JK      0.9448        0.0464      0.0002 
Weibull    down  LN      0.9433        0.0485      0.0002 


J. 


k  =  2 


Ditto 


k=1/2 


MP      0.9879         0.0691       0.0008 


JK      0.9222         0.1188      0.0087 


LN      0.8456        0.5116      0.0796 


MP      0.9727         0.5461       0.0854 


K.  Exponential  up  JK  0.9010  0.1032  0.0123 
Long-tailed  log  LN  0.9141  0.0817  0.0026 
normal   h    down        MP      0.9078        0.0792      0.0018 


L. 


h  =  0.2 


Ditto 


h=0.4 


JK       0.3835         0.1184       0.0183 


LN       0.9003         0.0918      0.0037 


MP      0.8930         0.0848       0.0024 


0.9454 
0.9300 
0.9422 

0.9394 
0.9228 
0.9336 

0.9489 
0. 9470 
0.9494 

0.9446 
0.7765 
0.8245 

0.9440 

0.9507 
0.9600 

0.9426 
0. 9511 
0.9604 
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Table   3:      Implementation  of   "Winsorizing"   on    Exponential   Up 
and   Samma    (k=1/2)     Down    Times    for    LNLJ    Procedure 

Sample    size=15 
Average  Variance  Point 

Coverage  Width  Width  Availability 


G=0 

0.9103 

0. 4212 

0.0809 

0.8239 

G=1 

0.9095 

0.3624 

0.0572 

0.8686 

G=2 

0.9272 

0.3467 

0.0567 

0.8840 

G=3 

0.9400 

0.3427 

0.0609 

0.8954 

Table   4:      Implementation  of   "Winsorizing"   on    Exponential   Up 
and   Gamma    (k=1/2)    Down   times    for    LNLJ    Procedure 

Sample   size=25 
Averaga  Variance  Point 

Coverage  Width  Width  Availability 


G=0 

0.8256 

0.3691 

0. 0670 

0.8319 

G=2 

0.3010 

0.2565 

0.0265 

0.8829 

G=3 

0.8361 

0.2456 

0.0247 

0.8908 

G=4 

0.8717 

0.2391 

0.0247 

0.8972 
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Table   5:      Implementation  of    "Winsorizing"   or.    Exponential   Up 
and   Weibull    (k=1/2)    Down   Times   for    LNLJ   Procedure 
Sample    size=15 

Average  Variance  Point 

Coverage  Width  Width  Availability 

G=0  0.9336  0.5752  0.0921  0.7649 

G=1  0.9418  0.5056  0.0828  0.8323 

G=2  0.9526  0.4699  0.0855  0.8611 

G=3  0.9510  0.4464  0.0907  0.8836 


Table  6:      Implementation  of   "Winsorizing"   on    Exponential   Up 
and   Weibull    (k=1/2)    Down   Times   for    LNLJ   Procedure 
Sample    size=25 
Average  Variance  Point 

Coverage  Width  Width  Availability 

G=0  0.8456  0.5116  0.0796  0.7765 

G=2  0.8692  0.3874  0.0516  0.8537 

G=3  0.9051  0.3632  0.0513  0.8689 

G=4  0.9316  0.3431  0.0511  0.8329 
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Table  7:      Implementation  of   "Biweighting"   on    Exponential   Op 

and   Gamma    (k=1/2)  Down    Times    for    LNLJ    Procedure 
Sample                          Average        Variance  Point 

size      Coverage  width  width  Availability 

Original  15  0.9103  0.4212  0.0809  0.8223 

Biweighted        15  0.8935  0.2621  0.0448  0.9117 

Original  25  0.8256  0.3691  0.0670  0.8319 

Biweighted        25  0.3786  0.1768  0.0159  0.9195 


Table   8:      Implementation   of   "Biweighting"   on  Exponential   Op 

and    Weibull    (k=1/2)     Down  Times    for  LNLJ   Procedure 

Sample                          Average  Variance           Point 

size      Coverage        width  width  Availability 

Original           15           0.9336           0.5752  0.0921  0.7649 

Biweighted         15           0.8423           0.3163  0.0784  0.9071 

Original           25           0.3456           0.5116  0.0796  0.7765 

Biweighted        25           0.8374           0.2134  0.0362  0.9231 
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Table   9:      Implementation  of   "Grouped  Jackknif e"    Procedure    for 

Exponential    Op   and   Weibull  (k=1/2)    Down   Times 

Sample                               Average  Variance               Point 

Method        Size        Coverage          Width  Width          Availability 


JK 


14 


0.9209 


3. 2136 


0.0353 


0.9381 


LN 


14 


0.9486 


0.3648 


0.0774 


0.9038 


MP 


14 


0.9390 


0.3083 


0.0672 


0.9073 


JK 


24 


0.9223 


0.  1319 


0.0118 


0.9429 


LN 


24 


0.9500 


0.2469 


0.0400 


0.9148 


MP 


24 


0.9685 


0.25  18 


0.0452 


0.9183 


Table    10:      Implementation   of  "Grouped    Jackknife"    Procedure 
for   Exponential    Up    and   Gamma    (k=1/2)    Down   Times 
Sample  Average        Variance  Point 

Method        Size        Coverage  Width  Width  Availability 


JK 

14 

0.9426 

0.1330 

0.0080 

0.9441 

LN 

14 

0.9534 

0.2266 

0.0381 

0.9243 

MP 

14 

0.9608 

0.2186 

0.0419 

0.9262 

JK 

24 

0.9472 

0.0839 

0.0015 

0.9461 

LN 

24 

0.9U62 

0  .  1  4  4  8 

0.0174 

0.9293 

MP 

24 

0.9783 

0.  1710 

0. 0259 

0.9322 
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Table  11:   Implementation  of  "Winsorizing"  on  Exponential  Up 

and  Long-Tailed  Log-Normal  (h=0.2)  Down  Times 
Sample  size=15 
Average     Variance       Point 
Coverage     Width       Width      Availability 


JK 

G=0 

0.8968 

0. 1385 

0.0220 

0.9416 

G=1 

0.8190 

0.0700 

0.0019 

0.9544 

G=2 

0.7450 

0.0585 

0.001 1 

0.9560 

G=3 

0.7160 

0.0  54 

0.0009 

0.9565 

LN 

G=0 

0.9119 

0.1217 

0.0091 

0.9483 

G=1 

0.8560 

0.0898 

0.0040 

0.9549 

G=2 

0.8040 

0.0759 

0.0029 

0.9575 

G=3 

0.7460 

0.0668 

0.0027 

0.9598 
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Table    12:      Implementation   of   "Winsorizing"    on  Exponential   Op 

and  Long-Tailed    Log-Normal     (h=0.2)  Down   Times 

Sample    size=25 
Average  Variance  Point 


Coverage 

Width 

Width 

Availabil 

JK 

G  =  0 

0.9010 

0.  1032 

0.0123 

0.9440 

G  =  2 

0.7380 

0.0460 

0.0004 

0.9554 

G=3 

0.6950 

0.0416 

0.0001 

0.9563 

G=4 

0.6620 

0.0337 

0.0003 

0.9568 

LN 

G  =  0 

0.9141 

0.0817 

0.0026 

0.9507 

G=2 

0.8240 

0.0594 

0.0010 

0.9558 

G  =  3 

0.7750 

0.0539 

0.0008 

0.9573 

G=4 

0.7380 

0.0483 

0.0006 

0.9583 
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Table    13:      Implementation   of   "Winsorizing"    on  Exponential   Up 

and  Long-Tailed   Log-Normal     (h=0.4)  Down   Times 
Sampls    size=15 
Average             Variance  Point 


JK 


Coverage 


Width 


Width 


Availability 


G=0 


0.8794 


0.1 578 


0.0305 


0.9394 


G=1 


0.7560 


0.0730 


0.0026 


0.9542 


G=2 


0  .6870 


0.0588 


0.0015 


0.9560 


G=3 


0.6460 


0.0537 


0.0012 


0.9565 


LH 


G=0 


G=1 


G=2 


G=3 


0.9009 


0.8470 


0.7800 


0.7180 


0. 1394 


0. 1049 


0.0866 


0.0  76  3 


0.0136 
0.0069 
0.0053 

0.0043 


0.9483 
0.9543 
0.9576 
0.9603 
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Table  14:   Implementation  of  '•» insorizing"  on  Exponential  Op 

and  Long-Tailed  Log-Normal  (h=0.4)  Down  Times 
Sample  size=25 
Average     Variance       Point 

Coverage     Width       Width  Availability 

JK 


G=0 

0.8835 

0.1  184 

0.0183 

0.9426 

G=2 

0.6810 

0.0473 

0.0007 

0.9560 

G=3 

0.6220 

0.0411 

0.0004 

0.9572 

G=4 

0.5750 

0.0380 

0.0004 

0.9578 

LN 

G=0 

0.9003 

0.0918 

0.0037 

0.951  1 

G=2 

0.8000 

0.0666 

0.0018 

0.9568 

G=3 

0.7660 

0.0585 

0.0014 

0.9590 

G=4 

0.7000 

0.0531 

0.0012 

0.9607 
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Table  15:   Implementation  of  "Biweighting"  on  Exponential  Up 
and  Long-Tailed  Log-Normal  (h=0.2)  Down  Times 
Sample  Average    Variance     Point 

size   Coverage   width     width    Availability 

JK 


Original 

15 

0.3968 

0.  1385 

0.0220 

0.9416 

Biweighted 

15 

0.6870 

0.0528 

0.0010 

0.9554 

Original 

25 

0.9010 

0. 1032 

0.0123 

0.9440 

Biweighted 

25 

0.6360 

0.0362 

0.0002 

0.9569 

LN 

Original 

15 

0.91  19 

0.  1217 

0.0091 

0.9483 

Biweighted 

15 

0.7010 

0.0604 

0.0019 

0.9607 

Original 

25 

0.9141 

0.0817 

0.0026 

0.9507 

Biweighted 

25 

0.6170 

0.0406 

0.0004 

0.9627 
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Table    16:      Implementation   of   "Biweighting"    on    Exponential   Up 
and  Long-Tailed   Log-Normal    (h=0.4)     Down   Times 
Sample  Average        Variance  Point 

size      Coverage        width  width  Availability 

JK 


Original 

15 

0.  87  94 

0.1578 

0.0305 

0.9394 

Biweighted 

15 

0.6210 

0.0U99 

0.3010 

0.9578 

Original 

25 

0.8835 

0.  1184 

0.0183 

0.9426 

Biweighted 

25 

0.5500 

0.0340 

0.0002 

0.9590 

LN 


Original 

15 

0.  90  09 

0. 1394 

0.0136 

0.9483 

Biweighted 

15 

0.6350 

0.0614 

0.0025 

0.9642 

Original 

25 

0.9003 

0.0918 

0.0037 

0.9511 

Biweighted 

25 

0.5261 

0.0394 

0.0005 

0.9663 
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IV.        CONCLUSIONS 

The    LNLJ    procedure    works    well    when    the    actual    down   times 
are   Exponential   or   Log-normal.      However,    the   procedure   is 
very   sensitive   to   the    sample   variance.      If    the   sample    vari- 
ance  tends    to  vary   greatly,    sometimes    becoming   excessively 
large,    then   the   method  tends   to   fail;    the    intervals   shifted 
down.      For    instance,    long    tailed   3amma,    weibull    and   Log-nor- 
mal  distributed   down   times    provide    extremely  large    sample 
variances,    and   under    these    circumtances   the   method    fails. 

For    further   study,    one    might   apply    jacxknifing    leaving 
out   one   observation   at   a   time,    instead   leaving   out    a    pair. 
Obviously,    this   is   going   to    increase   the    degrees   of    freedom, 
but   might   by   help   reducing    the    variance   and   the   confidence 
interval   length.      On   the   other    hand,    grouping   more    than   two, 
might   provide  more  stable    results. 

In    addition,    in    order    to   see   what    happens   in    one-sided 
availability   estimation,    Exponential    up   and   Log-normal   down 
and   Exponential   up   tfeibuil    (x=1/2)     down    times    were    simulated 
and   results   are   shown    in   Appendix   D. 
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APPENDIX  A 

COMPUTER  PROGRAMS 

Simulation  program  consists  of  main  program  and  several 
subroutines.   Main  program,  for  JK,  LNLJ  and  MP  procedures 
computes  availability  confidence  limits  and  point  availabil- 
ity and,  scores  the  coverage  for  each  replication.   Then, 
after  1000  replications  computes  the  statistics  of  these 
parameters  and  prints  out  the  results  for  given  underlying 
distributions. 

Subroutine  CONF  computes  two-sided  confidence  limits  of 
a  given  data  vector.   Subroutine  MOMENT  computes  the  fourth 
moment  of  given  data  vector.   Subroutine  LNLJ  generates  the 
up  and  down  times  from  given  underlying  distributions,  and 
takes  the  log  transform  of  these,  and  computes  the  means  and 
the  variances  of  these  data  vectors. 

Subroutine  BIWGT  computes  the  Biweight  estimates  of  a 
given  data  vector,  and  uses  Subroutine  MEDIAN  for  computing 
the  median  of  this  data.   Subroutine  WINSOR  Winsorize  a 
given  data  vector  for  given  "g'1  level. 
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APPENDIX    B 

IsQNG-IilLSD    LOG -NORMAL    DISTRIBUTION 

(1)  Let    W   be    log-normal    random   variable   which 
lnW/v^  N  ( ii,a2)  .      k   moment    of   W   as   follows, 
E[tfk]=   exp  (leu  +    1/2«k2    a2)  , 

see    [Ref.    9]. 

So, 

E[W]=    exp(   u  ♦    1/2*    a2) 

E[»2  ]=    exp  (2  m    2   a2) 

(2)  Create  a    stretched    log- normal    random   variable 
D=   c«W« (1+h«W) 

where  c   and   h   are   constant   and   0<h<1. 

E[D  ]=    c(E[  »]+h«E[w2  ]} 

E[D]=   c{exp(u+1/2«   a2)    ♦    h»exp  (2U    +2    a2)} 

c-    E[D]    {    exp  (  u  +1/2  •    a2)     +h «exp  (2  u  +2  a  2)  }  " l        (a.  1) 

(3)  On    the   other    andf    we   can   write    the    median    for    which 
is   the   function    of    the   expected   value   as    follows, 

M=    d  (0.5)  =  $  E[  D] 

M=    c«exp  (u)  •  (1+h«exp  (  \x )  )  =   *E[D]  (a.  2) 
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where  ^  is  the  median  location  factor  of  the  new  dis- 
tribution.  If  we  put  equation  (a.1)  in  equation 
(a.  2)  ,  we  obtain 

exp  (y)  •  (1  +  h«exp(u)) 

^•2[D]=  E[D] 

exp(y+1/2«    a2)     ♦    h«exp(2u  +  2    a2) 
After   the   cancellations,    we    obtain 

i|»  (exp(y  «■    1/2*  a2)     +h»exp  (2U   +    2    a2)    }=    1  +  h«exp  (y ) 
i\>  »exp  (1/2  a  2)    -    1    =    h«exp(u)  •  (1    -  ty  «exp  (2  a  2)  ) 
1-  <P  •exp(1/2«  a  2) 

h=   exp(-u>- — 

i|>   «exp(2  a  2)     -    1 
If   h  =  0,    then    ty   =  exp  (-1/2*   a2) 

1    +   exp  (u) 
if   h=1,    then     4>  = 


exp(u  +  2    a2)     ♦    exp  (1/2»    a2) 
So,    if  we    get 
exp(-1/2«a2)     <    ^<    exp(u*2    a2)  >exp  (1/2«    a2) 

then    we   obtain    h   between    0    and    1. 
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APPENDIX  S 

DISTRIBaTIONAL  FIGURES 
The  LNLJ  method  fails  for  Exponential  up  and  Gamma 
(k=1/2)  down,  Exponential  up  and  Weibull  (k=1/2)  down,  and 
Exponential  up  and  Long-tailed  Log-  normal  h  down  times 
cases.   In  order  to  demonsrate  the  distributional  charecter- 
istics  of  this  cases,  histograms  are  obtained  from  simula- 
tion for  up  times,  down  times,  log  transformations  of  down 
times  and  pseudovalues.   Also,  histograms  are  repeated  after 
grouping  the  observations  in  order  to  observe  the  grouping 
effect  on  the  parameters.   Distributional  charecteristics  of 
these  parameters  such  as  means,  variances,  moments  and 
shapes  can  be  observed  from  these  histograms. 
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APPENDIX    D 

ONE-SIDED    AVAILABILITY    ESTIMATION 
Best   and  worst   cases   for  LNLJ    procedure    were   simulated 
for   one-sided  availability    estimation,    and    results   are   shown 
in    Tables    17    and    18    as    follows; 


Table    17:      Exponential    Up   and   Log-normal   Down    Times 


Low.  Availability    Bound  Cover 

Upp.  Availability   Bound  Cover 

Avg.  Lower    Availability  3ound 

Avg.  Upper    Availability  Bound 


N=15 

JK  LN 


N=25 
JK  LN 


0.999       1.000      0.999       1.000 


0.933      0.899      0.966       0.965 


0.889       0.879      0.910       0.908 


0.966       0.964      0.966       0.966 


Table  18:   Exponential  Up  and  Weibull  (k=1/2)  Down  Times 


Low.  Availability  Bound  Cover 

Upp.  Availability  Bound  Cover 

Avg.  Lower  Availability  3ound 

Avg.  Upper  Availability  Bound 
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N=15 
JK  LN 


N=25 
JK  LN 


0.900       0.990      0.914       0.997 


0.966       0.853      0.955       0.705 


0.344       0.460      0.331        0.500 


0.972       0.970      0.975       0.955 


Notice    the    low    Average    Lower    Availability   Bound    provided 
by    the    LNLJ    procedure    for    the   tfeibull    (k=1/2).       Results   show 
that   the   one-sided   availability    estimation    behaves    as   the 
two-sided   availability   estimation. 
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